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Abstract. The dynamics of Gowdy vacuum spacetimes is considered in 
terms of Hubble-normalized scale-invariant variables, using the timelike area 
^ ■ temporal gauge. The resulting state space formulation provides for a simple 

| mechanism for the formation of "false" and "true spikes" in the approach 

£\j . to the singularity, and a geometrical formulation for the local attractor. 
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1. Introduction 

m ; 

Our map of general relativity has been, and still is, full of areas marked 
"Terra Incognita". Its oceans and forests are as full of monsters, as ever the 
■ outer limits of the world during the dark ages. Only recently hunters have gone 

CN . out to bring down some of the large prey of general relativity. This is a hunt 

that requires the utmost ingenuity and perseverance of the hunter. It is with 
pleasure that we dedicate this paper to Vince Moncrief, who has brought home 
a number of prize specimens, both in the field of general relativity as well as on 
the African veldt. 

The cosmic censorship conjecture is the most elusive of the game animals of 
general relativity. No one has come close to bringing it down, though it is the 
subject of many stories told around the camp fires. However, some of its lesser 
cousins, appearing as cosmologies with isometries, can and have been hunted 
successfully. Here, Vince Moncrief has consistently led the way. This is true for 
^ ■ spatially homogeneous (Bianchi) cosmologies, Gowdy vacuum spacetimes, as 

well as [/(l)-symmetric spacetimes. In addition, he has made fundamental con- 
tributions to our understanding of the structure of 3+1 dimensional spacetimes 
with no isometries. 

In this paper we will discuss Gowdy vacuum spacetimes, a class of spa- 
tially inhomogeneous cosmological models with no matter sources, two com- 
muting spacelike Killing vector fields, and compact spacelike sections. Gowdy 
spacetimes are named after Robert H. Gowdy who, inspired by a study of 
Einstein- Rosen gravitational wave metrics JB] \ initially investigated this class 
of metrics |22l 123) . The dynamics of Gowdy spacetimes exhibit considerable 
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complexity, and reflect the nonlinear interactions between the two polarization 
states of gravitational waves. The symmetry assumption restricts the spatial 
topology to S 3 ,S 1 xS 2 , or T 3 . In a series of papers, Moncrief and collaborators 
have built up a picture of the large scale structure of these spacetimes, and, in 
particular, of the nature of their cosmological singularities. 

In Ref. [SSI) Vince Moncrief laid the foundation for the analysis of Gowdy 
vacuum spacetimes, by providing a Hamiltonian framework and proving, away 
from the singularity, the global existence of smooth solutions to the associated 
initial value problem. The Hamiltonian formulation has played a central role in 
the later work with Berger; in particular, in the method of consistent potentials 
(see Refs. jl] and [H] for a discussion). 

Isenberg and Moncrief [2S] proved that for polarized Gowdy vacuum space- 
times (i) the approach to the singularity is "asymptotically velocity-term domi- 
nated" (AVTD), i.e., the asymptotic dynamics is Kasner-like, and (ii) in general 
the spacetime curvature becomes unbounded. Chrusciel et al constructed 
polarized Gowdy vacuum spacetimes with prescribed regularity at the singu- 
larity, while Isenberg et al [57j studied the evolution towards the singularity of 
the Bel-Robinson energy in the unpolarized case. 

Grubisic and Moncrief |24j studied perturbative expansions of solutions for 
Gowdy vacuum spacetimes with topology T 3 x R near the singularity, which are 
consistent when a certain "low velocity" assumption holds. Based on the intu- 
ition developed in this work, Kichenassamy and Rendall [2H] used the Fuchsian 
algorithm to construct whole families of solutions, depending on the maximum 
possible number of four (here analytic) free functions. 

Berger and Moncrief [E] performed numerical simulations of the approach 
to the singularity and found that "spiky features" form in the solutions near 
the singularity, and that these features "freeze in". A nongeneric condition 
used in the work of Grubisic and Moncrief |2l] turns out to be related to the 
formation of spikes. Further numerical experiments clarifying the nature of the 
spiky features have been carried out by Berger and Garfinkle 0, and recently 
by Garfinkle and Weaver [2^. Garfinkle and Weaver studied the dynamical 
behavior of general spikes with "high velocity" in Gowdy vacuum spacetimes 
with topology T 3 x R, using a combination of a Cauchy and a characteristic 
integration code. They found numerical evidence of a mechanism that drives 
spikes with initially "high velocity" to "low velocity" so that only spikes of the 
latter type appear to persist in the approach to the singularity. 

The understanding of the spiky features for generic Gowdy vacuum space- 
times is the main obstacle to a complete understanding of the nature of the 
singularity in this setting. Rendall and Weaver [31] have constructed families 
of Gowdy vacuum spacetimes with spikes by an explicit transformation, starting 
from smooth "low velocity" solutions. In addition, they introduced a classifica- 
tion of spikes into "false spikes" (also: "downward pointing spikes") and "true 
spikes" (also: "upward pointing spikes"), which by now has become part of the 
standard terminology in this area. Substantial progress in understanding the 
structure of Gowdy vacuum spacetimes going beyond what is possible by means 
of the Fuchsian algorithm has been made recently by Chae and Chrusciel [0] 
and Ringstrom [SHI EZ[ • 



GOWDY PHENOMENOLOGY IN SCALE-INVARIANT VARIABLES 



3 



In spite of the fact that Hamiltonian methods are Vince's weapon of choice, 
in this paper we will bring the dynamical systems formalism, developed by 
John Wainwright, Claes Uggla and collaborators, to bear on the problem of 
understanding the dynamics of Gowdy vacuum spacetimes near the singular- 
ity (assuming spatial topology T 3 throughout). This approach is a comple- 
mentary alternative to the method of consistent potentials, often used in the 
work of Moncrief, Berger and collaborators. The dynamical systems formalism 
breaks the canonical nature of the evolution equations by passing to a system 
of equations for typically Hubble-normalized scale-invariant orthonormal frame 
variables. The picture of a dynamical state space with a hierarchical structure 
emerges naturally. By casting the evolution equations into a first order form, 
this approach allows one to extract an asymptotic dynamical system which ap- 
proximates the dynamics near the singularity increasingly accurately. The guid- 
ing principle for understanding the approach to the singularity is "asymptotic 
silence" , the gradual freezing-in of the propagation of geometrical information 
due to rapidly increasing spacetime curvature. One finds that in this regime, 
with spacetime curvature radii of the order of the Planck scale, the dynamics 
along individual timelines can be approximated by the asymptotic dynamics 
of spatially self-similar and spatially homogeneous models — a feature directly 
associated with the approach to the so-called "silent boundary." 

1.1. Dynamical systems methods. In this paper, we will study Gowdy 
vacuum spacetimes in terms of Hubble-normalized scale-invariant orthonormal 
frame variables. This approach was first pursued systematically in the context 
of spatially homogeneous (Bianchi) cosmology by Wainwright and Hsu |3J 2 , 
building on earlier work by Collins ^2]- The main reference on this material is 
the book edited by Wainwright and Ellis (WE) j^DJ , where techniques from the 
theory of dynamical systems are applied to study the asymptotic dynamics of 
cosmological models with a perfect fluid matter source that have isometries. 

Hewitt and Wainwright [201 obtained a Hubble-normalized system of equa- 
tions for spatially inhomogeneous cosmological models with a perfect fluid mat- 
ter source that admit an Abelian G2 isometry group which acts orthogonally 
transitively (i.e., there exists a family of timelike 2-surfaces orthogonal to the 
GVorbits). 3 When restricting these models to the vacuum subcase, imposing 
topolog y T 2 on the GVorbits, and compactifying the coordinate associated with 
the spatial inhomogeneity (so that overall spatial topology T 3 results), one gets 
Gowdy vacuum spacetimes as an invariant set. 

The nature of the singularity in orthogonally transitive perfect fluid G2 cos- 
mologies has been studied by dynamical systems methods in Ref. [20], where 
the area expansion rate of the G^-orbits was used as a normalization variable 
instead of the more common Hubble scalar H. In recent work, Uggla et al [3EI 
introduced a general dynamical systems framework for dealing with perfect 
fluid Go cosmologies, i.e., the full 3+1 dimensional case without isometries. 

2 More accurately, in this work the authors employed the volume expansion rate O as the 
normalization variable rather than the Hubble scalar H; these two variables are related by 
9 = 3H. 

3 See previous footnote. 
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Using again a Hubble-normalized system of equations, this work studies the as- 
ymptotic dynamics of the approach to the singularity. Physical considerations 
suggest that the past directed evolution gradually approaches a boundary of 
the Go state space on which the Hubble-normalized frame variables E a l , i.e., 
the coefficients of the spatial derivatives, vanish. This so-called "silent bound- 
ary" is distinguished by a reduction of the evolution equations to an effective 
system of ordinary differential equations, with position dependent coefficients. 
In this dynamical regime, typical orbits in the Go state space shadow more and 
more closely the subset of orbits in the silent boundary, governed by the scale- 
invariant evolution equations restricted to the silent boundary, which coincide 
with the evolution equations for spatially self-similar and spatially homogeneous 
models. It is a strength of the dynamical systems framework that it provides 
for the possibilities (i) to make some interesting conjectures on the nature of the 
local past attractor for the Hubble-normalized system of equations of Go cos- 
mology, and, in particular, (ii) to give a precise mathematical formulation of the 
well known BKL 4 conjecture in theoretical cosmology: For almost all cosmolog- 
ical solutions of Einstein's field equations, a spacelike initial singularity is silent, 
vacuum dominated and oscillatory (see Refs. §§2,3], [HI §§3,5] and [3E1 §6]). 
Both the AVTD behavior of Gowdy vacuum spacetimes and the past attractor 
associated with the BKL conjecture turn out to be associated with invariant 
subsets on the silent boundary. Moreover, the dynamics on the silent boundary 
also governs the approach towards these subsets. Hence, the silent boundary 
provides a natural setting for studying AVTD as well as asymptotic oscillatory 
BKL behavior, and how such behavior arises. 

1.2. Overview of this paper. In Sec.[5J we introduce the system of equations 
for Gowdy vacuum spacetimes with topology T 3 x R in the metric form used 
in the work of Moncrief and collaborators, and also describe the orthonormal 
frame variables. In Sec. 01 we introduce two alternative scale-invariant systems 
of equations. In Subsec. 13. 1[ we introduce Hubble-normalized variables, and 
derive the system of evolution equations and constraints which govern the dy- 
namics of the Hubble-normalized state vector X(t, x). In Subsec. 13.21 on the 
other hand, we discuss the area expansion rate-normalized variables and equa- 
tions, introduced in Ref. [201. Invariant sets of the Hubble-normalized system 
of equations for Gowdy vacuum spacetimes are described in Sec. 0J the most 
notable ones being the polarized models and the silent boundary of the Hubble- 
normalized state space; the latter is associated with spatially self-similar and 
spatially homogeneous models (including the Kasner subset). In Sec.JSl we start 
the discussion of the asymptotic dynamics towards the singularity of the Gowdy 
equations. In Subsec. [5. 1\ we review the picture found numerically by Moncrief, 
Berger and collaborators, and in Subsec. ff>.2l we describe the phenomenology in 
terms of Hubble-normalized variables, using principles introduced in Refs. jSH] 
and |2U| . In particular, we numerically show that there is good agreement be- 
tween state space orbits of the full Gowdy equations for X(t,x) and state space 
orbits of the system of equations on the silent boundary. We thus give evidence 
for the hypothesis that the silent boundary describes both the final state and 
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the approach towards it, as claimed above. We then derive and solve the sys- 
tem of ordinary differential equations which governs the asymptotic behavior 
of the state vector X(t, x) in a neighborhood of the Kasner subset on the silent 
boundary. The connection between these analytic results and others obtained 
through various numerical experiments are subsequently discussed. In partic- 
ular, we point out the close correspondence between the formation of spikes 
and state space orbits of X(t,x) of Bianchi Type-I and Type-II, for timelines 
x = constant close to spike point timelines. In Subsec. 15.31 we then comment 
on the behavior of the Hubble-normalized Weyl curvature in the "low velocity" 
regime. In Sec. we conclude with some general remarks. We discuss the 
dynamical nature of the silent boundary of the state space, relate this concept 
to the AVTD and BKL pictures in Subsec. 16.11 and review its application to 
Gowdy vacuum spacetimes in Subsec. 16.21 We also comment on the conjec- 
tural picture that emerges from the present work, and make links to a more 
general context. Finally, expressions for the Hubble-normalized components of 
the Weyl curvature tensor, as well as some scalar quantities constructed from 
it, are given in an appendix. 



2.1. Metric approach. Gowdy vacuum spacetimes with topology T 3 x R are 



using the sign conventions of Berger and Garfinkle [Hj, but denoting by t their 
time coordinate r. The metric variables A, P and Q are functions of the local 
coordinates t and x only. Each of the spatial coordinates x, y\ and 1/2 has 
period 2ir, while the (logarithmic) area time coordinate t runs from — 00 to 
+00, with the singularity at t = +00. This metric is invariant under the 
transformations generated by an Abelian Gi isometry group, with spacelike 
Killing vector fields £ = d yi and rj = d y2 acting orthogonally transitively on T 2 
(see, e.g., Hewitt and Wainwright [2E])- Note that the volume element on T 3 
is given by y^g = £qC _ ( 3 * +a )/ 4 . 




Einstein's field equations in vacuum give for P and Q the coupled semilinear 
wave equations 



2. Gowdy vacuum spacetimes 





P,tt ~ e~ 2t P, xx = e 2P (Q 2 - e~ 2t Q 2 x ) , 
Q,tt - e~ 2t Q, xx = - 2 (P t Q >t - e~ 2t P x Q tX ) 



(2.2a) 
(2.2b) 



and for A the Gaufi and Codacci constraints 



= A, - [P 2 + e~ 2t P 2 x + e 2P {Q 2 t + e~ 2t Q 2 x )] , 



(2.3a) 
(2.3b) 




^The most general form of the Gowdy metric on T 3 x R has been given in Ref. |1UI §2]. 
However, the additional constants appearing in that paper have no effect on the evolution 
equations. 
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For Gowdy vacuum spacetimes the evolution (wave) equations decouple from 
the constraints. Consequently, the initial data for P, Q, Pi and as 2ir- 
periodic real-valued functions of x can be specified freely, subject only to ap- 
propriate requirements of minimal differentiability and the zero total linear 
momentum condition [33] 

= / X )X dx . (2.4) 
Jo 

The fact that the initial data contains four arbitrary functions supports the in- 
terpretation of Gowdy vacuum spacetimes as describing nonlinear interactions 
between the two polarization states of gravitational waves. The so-called polar- 
ized Gowdy vacuum spacetimes form the subset of solutions with = Q,t = Q, x - 
Note that in this case, without loss of generality, a coordinate transformation on 
H 2 can be used to set = Q identically, and thus obtain the standard diagonal 
form of the metric. 

As is well known, the metric functions P and Q can be viewed as coordinates 
on the hyperbolic plane H 2 (this is just the Teichmiiller space of the flat metrics 
on the T 2 orbits of the G2 isometry group), with metric 

<? H 2 = dP 2 + e 2P dQ 2 . (2.5) 

The Gowdy evolution equations form a system of wave equations which is essen- 
tially a 1+1 dimensional wave map system with a friction term, the target space 
being H 2 . The sign of the friction term is such that the standard wave map 
energy decreases in the direction away from the singularity, i.e., as t — > —00. 
The friction term diverges as t — > +00. This is more easily seen using the 
dimensional time coordinate r = ^o e_t - I n terms of this time coordinate, the 
wave equations take the form 

-d 2 T u--d T u + dlu + T(du-du) = , (2.6) 

T 

where this can be seen explicitly. Here u is a map u: R 1,1 — > H 2 , and T are 
the Christoffel symbols defined with respect to the coordinates used on H 2 . It 
is convenient to define a hyperbolic velocity v by 

v:=\\d t u\\ u , = ^/p 2 + e^Q 2 t . (2.7) 

Then v(t,x) is the velocity in (H 2 ,g H 2) of the point [P(t,x), Q(t,x)] (cf. 
Ref. 

2.2. Orthonormal frame approach. In the orthonormal frame approach, 
Gowdy vacuum spacetimes arise as a subcase of vacuum spacetimes which admit 
an Abelian G2 isometry group that acts orthogonally transitively on spacelike 2- 
surfaces (cf. Refs. |26| and HO]). One introduces a group invariant orthonormal 
frame { e a } a =o,i,2,3) with eo a vorticity-free timelike reference congruence, e\ 
a hypersurface-orthogonal vector field defining the direction of spatial inhomo- 
geneity, and with e.2 and e% tangent to the GVorbits; the latter are generated by 
commuting spacelike Killing vector fields £ and r\ (see Ref. |391 §3] for further 
details). It is a convenient choice of spatial gauge to globally align e.2 with 
Thus, within a (3+l)-splitting picture that employs a set of symmetry adapted 
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local coordinates x M = (t,x,yi,y2), the frame vector fields e a can be expressed 
by 

e = N' 1 d t , e 1 = &i d x , e 2 = e 2 2 d yi , e 3 = e 3 2 d yi + e 3 3 d y2 , (2.8) 

where, for simplicity, the shift vector field N l was set to zero. 

For Gowdy vacuum spacetimes with topology T 3 x R, the orthonormal frame 
equations (16), (17) and (33)-(44) given in Ref. ^5] need to be specialized by 
setting (a, = 1, 2, 3) 

= e 2 (/) = e 3 (/) , (2.9a) 

= uj a = cr 3 i = a 12 = a a = n la = n 33 = ii 2 = u 3 = f2 2 = ^3 , (2.9b) 

= = p = q a = -K aP , (2.9c) 

for / any spacetime scalar. The remaining nonzero connection variables are 

(1) 0, which measures the volume expansion rate of the integral curves 
of eo; it determines the Hubble scalar H by H := | 0, 

(2) (Tn, cr 22 , <7 33 and cr 23 (with = an + cr 22 + cr 33 ), which measure the 
shear rate of the integral curves of eo, 

(3) n22 and n 23 , which are commutation functions for e±, e 2 and e 3 and 
constitute the spatial connection on T 3 , 

(4) iii, the acceleration (or nongeodesity) of the integral curves of eo, and 

(5) Qi, the angular velocity at which the frame vector fields e 2 and e 3 rotate 
about the Fermi propagated 6 frame vector field e±. 

The area density A of the G 2 -orbits is defined (up to a constant factor) by 

A 2 := (U a )(VbV b ) ~ (W 1 ) 2 , (2.10) 

which, in terms of the coordinate components of the frame vector fields e 2 and 
e 3 tangent to the G 2 -orbits, becomes 

A- 1 = e 2 2 e 3 3 . (2.11) 

The physical dimension of A is [length] 2 . The key equations for A, derivable 
from the commutators (16) and (17) in Ref. |19j . are 

N- 1 ^ = (2H - an) , ei^ = 0, (2.12) 

i.e., (2H — an) is the area expansion rate of the G 2 -orbits. 

We now introduce connection variables a + , <r„, a x , Ti- and n x by 

a + = 2 ( U22 + = ~ 2 ' a - = 2~y| ( CJ 22 - 0-33) , (2.13a) 

1 1 

= -j= &23 , n- = — y= n 2 2 , (2.13b) 

1 

n x = -j= n 23 . (2.13c) 

Together with the Hubble scalar H, for Gowdy vacuum spacetimes this is a 
complete set of connection variables. 



spatial frame vector field e a is said to be Fermi propagated along eo if (e^, Ve e„) = 0, 
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Recall that e2 was chosen such that it is globally aligned with One finds 
that Qi = —\/3a x (see Ref. [201 §3]). The spatial frame gauge variable Q\ 
quantifies the extent to which e.2 and e 3 are not Fermi propagated along eo 
in the dynamics of Gowdy vacuum spacetimes, starting from a given initial 
configuration. 

We now give the explicit form of the natural orthonormal frame associ- 
ated with the metric (|2.1|) . With respect to the coordinate frame employed 
in Eq. (|2.1j) . the frame vector fields are 



eo 


= n- 1 d t 


= Ve (3t+A)/4 5 t , 


(2.14a) 


ei 


= ei 1 d x 


= l^e-«-W d x , 


(2.14b) 


e-2 


= £2 dyi 




(2.14c) 


e 3 


= e 3 2 d yi + e 3 3 d y2 


= i^e^/ 2 (-Qdy, + dy 2 ) . 


(2.14d) 



Note that the lapse function is related to the volume element on T 3 by N = 
£q 2 \/^g- Hence, in addition to being an area time coordinate, t is simultaneously 
a wavelike time coordinate. 

3. Scale-invariant variables and equations 

3.1. Hubble-normalization. Following WE and Ref. [3*5] . we now employ the 
Hubble scalar H as normalization variable to derive scale-invariant variables 
and equations to formulate the dynamics of Gowdy vacuum spacetimes with 
topology T 3 x R. Hubble-normalized frame and connection variables are thus 
defined by 

(Af-\E 1 1 ):=(N-\e 1 1 )/H (3.1) 
(£..., N.., U, R) := (<T..., n..., u x , Sli)/H . (3.2) 

In order to write the dimensional equations 8 in dimensionless Hubble-normalized 
form, it is advantageous to introduce the deceleration parameter, q (which is 
a standard scale- invariant quantity in observational cosmology), and the loga- 



rithmic spatial Hubble gradient, r, by 

(g + 1) : = -Af^dtHioH) , (3.3) 
r := - ft 1 d x \n(£oH) . (3.4) 
These two relations need to satisfy the integrability condition 

AT 1 d t r - 8 x q =(q + 2S+) r - (r - U) {q + 1) , (3.5) 

By definition, wavelike time coordinates satisfy the equation V M Vi/t = /(a; M ); / is 



a freely prescribable coordinate gauge source function of physical dimension [length] -2 . In 
terms of the lapse function, this equation reads N^ 1 (dt — N' di)N = 3HN + fN 2 . For the 
metric (12. 11 we have N 1 = and / = 0; the lapse equation thus integrates to oc In 
the literature, wavelike time coordinates are often referred to as "harmonic time coordinates" . 

8 In the orthonormal frame formalism, using units such that G/c 2 = 1 = c, the dynami- 
cal relations provided by the commutator equations, Einstein's field equations and Jacobi's 
identities all carry physical dimension [length] -2 . 
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as follows from the commutator equations (see Ref. [3E| for the analogous 
Go case). The Hubble- normalized versions of the key equations (|2.12j) are 

A/- 1 ^ = 2(l + E + ), £1^ = 0, (3.6) 

so that the magnitude of the spacetime gradient V a A is 

{V a A){V a A) = -AH 2 {l + Y, + fA 2 . (3.7) 

Thus, V a A is timelike as long as (1 + E + ) 2 > 0. As discussed later, £ + = — 1 
yields the (locally) Minkowski solution. 

In terms of the Hubble-normalized variables, the hyperbolic velocity, defined 
in Eq. (|2~7|) . reads 



hence, it constitutes a scale-invariant measure of the transverse (with respect 
to e±) shear rate of the integral curves of eg. 

3.1.1. Gauge choices. The natural choice of temporal gauge is the timelike area 
gauge (see Ref. |2U1 §3] for terminology). In terms of the Hubble-normalized 
variables, this is given by setting 

= + C , CGR\{0}, (3.9) 

N l = 0, so that by Eqs. (|3.6|) we have 

A = i\e 2Ct , (3.10) 

implying eo || V «4. This choice for TV and N l thus defines the dimensionless 
(logarithmic) area time coordinate t. For convenience we keep the real-valued 
constant C, so that the direction of t can be reversed when this is preferred. 
With the choice C = — ^, t coincides with that used in the metric (|2.1|) . 

Having adopted a spatial gauge that globally aligns e2 with £ has the conse- 
quence (see Ref. j^Ul §3]) 

R = -\/3S x . (3.11) 

This fixes the value of the Hubble-normalized angular velocity at which €.2 
and 63 rotate about e.\ during the evolution. Thus, while in this spatial gauge 
ei is Fermi propagated along eo, e.2 and e% are not. 

3.1.2. Hubble-normalized equations. The autonomous Hubble-normalized con- 
straints are 

(r — U) = -Ei 1 d x ln(l + £+) , (3.12a) 

1 = £;]_ + E 2 . + £ 2 X + iV^ +N 2 , (3.12b) 

(l + E + )?7 = -3(7V X £_ -iV_E x ) . (3.12c) 

These follow, respectively, from the commutators and the Gaufi and the Codacci 
constraints. The first is a gauge constraint, induced by fixing M according 
to Eq. (|in3>- Note that Eq. (I3.12bl) bounds the magnitudes of the Hubble- 
normalized variables E + , E_, E x , N x and N_ between the values and 1. 



10 L. ANDERSSON, H. VAN ELST, AND C. UGGLA 

The autonomous Hubble-normalized evolution equations are 

C~\l + £+) OtE^ = (g + 2£+) Ei 1 , (3.13a) 

C~\l + £+) $(1 + £+) = (q - 2) (1 + £+) , (3.13b) 
C-^l + £+) d t £_ + T^ 1 d x N x = (q - 2) £_ + (r - 17) iV x 

+ 2\/3£x - 2\/3 AT^ , (3.13c) 
C-^l + £+) d t N x + E x l c\£_ = {q + 2£ + ) N x + (r - 17) £_ , (3.13d) 
C- X (l + £+) 3 t S x - E 1 ! 1 fi^JV- = (g - 2 - 2^3£-) £ x 

- (r - U + 2\/3 AT X ) N- , (3. 13e) 
C _1 (l + £+) d t AT_ - T^ 1 9 X E X = (g + 2£ + + 2^3£_) AT_ 

-{r-U - 2v / 3A r x ) S x , (3.13f) 

where 

g = 2(£ 2 + E 2 . + £ 2 X ) - i (Ex 1 d x — r + U)U . (3.14) 

Note that from Eq. (I3.13bl) g can also be expressed by 

g = 2 + C- 1 ^(l + S+) . (3.15) 

Equations (|3.12|) - (|3.14j) govern the dynamics on the Hubble-normalized state 
space for Gowdy vacuum spacetimes of the state vector 

X = (# l 1 ,£ + ,£_,JV x ,£ x ,jV_) T . (3.16) 

Initial data can be set freely, as suitably differentiable 2-7r-periodic real-valued 
functions of x, for the four variables £_, iV x , £ x and iV_. Up to a sign, the 
data for £ + then follows from the Gaufi constraint l)3.12b|) . while the data 
for U follows from the Codacci constraint (|3,12c|) . The data for r, finally, is 
obtained from the gauge constraint (|3.12a[) . Note that there is no free function 
associated with E\ 1 . Any such function can be absorbed in a reparametrization 
of the coordinate x. 

It is generally helpful to think of the variable pair (£_,iV x ) as relating to 
the "+- polarization state" and the variable pair (£ x ,iV_) as relating to the 
"x -polarization state" of the gravitational waves in Gowdy vacuum spacetimes 
(see also Ref. [SB p. 58]). 

3.1.3. Conversion formulae. We will now express the Hubble-normalized frame 
and connection variables in terms of the metric variables A, P and Q. For the 
metric (|2.1|) . the choice for the nonzero real- valued constant C is 

The Hubble scalar, given by H = | N~ 1 (dt\/^g)/ y/^g, evaluates to 

= _£-l e (3i+A)/4 3+_\t (317) 
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As, by Eq. (|2,3a|) . A ;t > 0, we thus have H < (i.e., the spacetime is contracting 
in the direction t — ► +oo). We then find 

1 1 



(AT 1 , E^) 



U 



A x ) 



(E x , JV-) 



12 

3- 



3 + A, t 



3 + A, 



3 + A )t 
2V3 



3e~ 
Pi 



A 



3 + A + 



P, 



3 + A.^ 



2^3 [ 



3 + A, t 
-(t-P) 



Q,. 



V 3 + A,; 



3 + A, ; 



(3.18a) 
(3.18b) 
(3.18c) 
(3.18d) 



and R = - \/3E x . Note that Eq. (|3.18a|) leads to AfEi 1 = e~*, implying that 
for the metric (|2.1[) the timelike area and the null cone gauge conditions are 
simultaneously satisfied (see Ref. |2U1 §3] for terminology). The latter corre- 
sponds to choosing t to be a wavelike time coordinate. 

The auxiliary variables q and r, defined in Eqs. (|3.3|) and (|3.4[) . are given by 

12A,« _ t 3 / 4A, ix 



2 + 



3 + A* 



A.- + 



3 + A. ; 



(3.19) 



(3 + A,t) 2 ' 

Note that the gauge constraint (j3.12a() reduces to an identity when r, U, E\ l 
and E + are substituted by their expressions in terms of A, P and Q, given in 
Eqs. (|3. 19|) and (|3.18|) , Thus, this gauge constraint is redundant when formulat- 
ing the dynamics of Gowdy vacuum spacetimes in terms of the metric variables 
of Eq. (|2.1|) : it has been solved already. 

3.2. Area expansion rate normalization. Noting that the variable combi- 
nation (r — U), as well as the deceleration parameter q, show up frequently on 
the right hand sides of the Hubble- normalized evolution equations (|3.13|) . and 
so making explicit use of the gauge constraint (|3.12aj) and Eq. 1)3.15(1 . one finds 
that (for E + 7^ — 1) it is convenient to introduce variables E\ , E_, E X) A_ 
and A x by 

(£i\E_,Ex,A_,Ax) := (£ 1 1 ,E_,E X ,A_,A X )/(1 + E+) . (3.20) 

In terms of these variables, Eqs. (|3.13j) convert into the unconstrained au- 
tonomous symmetric hyperbolic evolution system 

C~ l d t E x l = 2E x l , 
C -1 $E_ + Ex 1 d x N x = 2Vst 2 x - 2V3N* 



C~ 1 d t N x +E l 1 d x t^ = 2N x , 

C~ l d t t x - Ek 1 d x N- = - 2v / 3E„ E x - 2y/3N x A_ , 
C~ l d t N- - Ei l d x t x = 2A_ + 2v / 3S_ A_ + 2y/3N x E> 
with an auxiliary equation given by 



c~ l d t 



1 



(1 + ^ 



3 



u 



(l + E. 



2N% + 2Nl 



(3.21a) 
(3.21b) 
(3.21c) 
(3.21d) 
(3.21e) 

(3.22) 
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The latter is redundant though, since the variables U := 17/(1 + £+) and 
S + := S+/(l + £+) are given algebraically by 

17 = f , (3.23a) 

S+ = i (1 - £ 2 - iV 2 - £ 2 - TV 2 ) , (3.23b) 

r = -3(iV x £_ -iV_£ x ) , (3.23c) 

the relations to which Eqs. 1)3.12(1 transform. But Eqs. (|3,21|) - (|3.23j) are nothing 
but the Gowdy equations written in terms of the area expansion rate normalized 
variables of Ref. [20]; it is obtained by appropriately specializing Eqs. (113)— 
(117) in that paper. The evolution system is well defined everywhere except (i) 
at X+ = — 1, and (ii) (for C > 0) in the limit t — > +oo. 

We now give the area expansion rate normalized frame and connection vari- 
ables in terms of the metric variables A, P and Q, and C = — \. Thus 

AT 1 , eA = - 2 (1, e -t ) , (3.24a) 



t+, V) =Ul-l\ tU e -*xA , (3.24b) 

7V X ) = (- P, t , e- l P x ) , (3.24c) 

(S x , 7V_) = - (e p Q, t , e"^ Q jS ) , (3.24d) 

and ^ = - v^Sx (see also Ref. |23 §A.3]). 

We will employ the Hubble-normalized variables and equations in the discus- 
sion which follows. 



4. Invariant sets 

For the Hubble-normalized orthonormal frame representation of the system 
of equations describing the dynamics of Gowdy vacuum spacetimes, given by 
Eqs. (|3.12|) - (|3.14j) . it is straightforward to identify several invariant sets. The 
most notable ones are: 

(1) Gowdy vacuum spacetimes with 

= S x = iV_ , 

referred to as the polarized invariant set, 

(2) Vacuum spacetimes that are spatially self-similar (SSS) in the sense of 
Eardley ^3 §3], given by 

= 8 X X , (7 = r , 

and 17 and r as in Eqs. (|3.12c|) and (|3.4[) . We refer to these as the SSS 
invariant set. Within the SSS invariant set we find as another subset 
spatially homogeneous (SH) vacuum spacetimes, obtained by setting 
= U = r. These form the SH invariant set. 
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(3) the invariant set given by 

= E x x , 

which, following Ref. |381 §3.2.3], we refer to as the silent boundary. 

Note that for the plane symmetric invariant set, defined by = £ x = iV_ = 
E_ = iVx , one has £+ = ±1 due to the Gaufi constraint (|3. 12bf) : E + = —1 yields 
a locally Minkowski solution, while E + = 1 yields a non-flat plane symmetric 
locally Kasner solution. Thus, contrary to some statements in the literature, 
general Gowdy vacuum spacetimes are not plane symmetric. 

There exists a deep connection between the SSS invariant set and the silent 
boundary: in the SSS case the evolution equation 1)3.1 3a|) for E\ l decouples 
and leaves a reduced system of equations for the remaining variables 9 ; the 
evolution equations in this reduced system are identical to those governing the 
dynamics on the silent boundary, obtained by setting Ei 1 = 0. The reduced 
system thus describes the essential dynamical features of the SSS problem. In 
contrast to the true SSS and SH cases, the variables on the silent boundary are 
x-dependent; thus orbits on the silent boundary do not in general correspond 
to exact solutions of Einstein's field equations. As will be shown, the silent 
boundary plays a key role in the asymptotic analysis towards the singularity 
of the Hubble-normalized system of equations. Therefore, we now take a more 
detailed look at the silent boundary and some of its associated invariant subsets. 

4.1. Silent boundary. Setting Ei 1 = in the Hubble-normalized system of 
equations (|3.12|) - (|3.14|) describing Gowdy vacuum spacetimes (and hence can- 
celling all terms involving E\ l d x ), reduces it to the system of constraints 

( r - U) = , (4.1a) 

1 = Y? + + E 2 . + E 2 X +Nl + iV 2 , (4.1b) 

(l + £ + )[/ = -3(7V x £_-7V_£ x ) , (4.1c) 

and evolution equations 

C-\l + E + ) d t {l + £+) = (q - 2) (1 + £+) , (4.2a) 

C _1 (l + E+) <%£_ = (q - 2) E_ + 2v / 3E 2 x - 2V3N 2 . , (4.2b) 

C-\l + £+) 8 t N x = (q + 2E+) N x , (4.2c) 

C^il + £+) d t ^ x = (q-2- 2v / 3S_) E x - 2y/3 N x N_ , (4.2d) 

CT\l + E+) d t N- = (q + 2E+ + 2a/3E_) iV_ + 2y/3N x E x , (4.2e) 

with 

g = 2(E 2 ^ + E 2 _ + E 2 x ) . (4.3) 
The integrability condition (|3.5|) reduces to 

C-\l + E + )d t r= (g + 2E + )r . (4.4) 



^The variable Ei 1 only appears as part of the spatial frame derivative E\ 1 d x in the equa- 
tions for these variables; the SSS requirement implies that these terms have to be dropped. 
We discuss these issues again in Sec.|Sl 
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In the above system, the variables U and r decouple, and so Eqs. (|4.1a|) . (|4.1c|) 
and Q4.4J) can be treated separately; the remaining autonomous system thus 
contains the key dynamical information. Note that it follows from Eqs. (|4.2cj) 
and (|4.4fl that r is proportional to N x , with an x-dependent proportionality 
factor. This is not a coincidence: for the corresponding SSS equations the 
square of this factor is proportional to the constant self-similarity parameter /, 
defined in Ref. pH p. 294]. 

On the silent boundary we see a great simplification: firstly, the system 
is reduced from a system of partial differential equations to one of ordinary 
differential equations; secondly, the reduction of the gauge constraint l)3.12ajl 
to Eq. (|4.1aj) yields a decoupling of U and r; note, however, that the Gaufi and 
Codacci constraints, Eqs. (|3.12b|) and (|3.12c|) . remain unchanged. 

We will see in Sec.EJbelow that the "silent equations," Eqs. (|3.12b|) . l|3,12cj) . 
(|4.2|) and (|4.3j) . significantly influence and govern the asymptotic dynamics of 
Gowdy vacuum spacetimes in the approach to the singularity. 

4.1.1. Polarized set on the silent boundary. Taking the intersection between 
the polarized invariant set (0 = S x = iV_) and the silent boundary (0 = -Ei 1 ) 
yields: 

1 = T? + + + N 2 X , (4.5a) 
(1 + £+) U = - 3S„ iV x = (l + S+)r, (4.5b) 



and 



C- 1 d t (l + £+) = -2Nl , (4.6a) 



C-\l + Y, + )d t N x =2 



N x , (4.6b) 



£_ = constant x (1 + E + ) . (4.6c) 

Note that with C < 0, N x increases to the future when (S + ,S_) take val- 
ues inside a disc of radius \ centered on (—^,0) in the (S + S_)-plane of the 
Hubble-normalized state space, and decreases outside of it. Note also that the 
projections into the (S + S„)-plane of the orbits determined by these equations 
constitute a 1-parameter family of straight lines that emanate from the point 
(—1,0) (which we will later refer to as the Taub point T\). 

4.2. SH invariant set. By setting = U = r one obtains the SH invariant set 
as a subset of the SSS invariant set. Since in the present case n a p n a/3 — (n a a ) 2 > 
0, it follows that this invariant set contains vacuum Bianchi models of Type-I, 
Type-II and Type-VIo only (see King and Ellis |3U1 p. 216]), though in general 
not in the canonical frame representation of Ellis and MacCallum |17j (presently 
the spatial frame { e a } is not an eigenframe of N a p, nor is it Fermi propagated 
along eo). However, the conversion to the latter can be achieved by a time 
dependent global rotation of the spatial frame about e%. Note that the present 
representation means that there exists only one reduction (Lie contraction) from 
Type-VIo to Type-II, in contrast to the diagonal representation where there 
exist two possible reductions. This has the consequence that the dynamics 
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of Gowdy vacuum spacetimes allows only for the existence of a single kind of 
curvature transition (see Sec.EJ), described by the Taub invariant set introduced 
below. 

The SH restrictions, = d x X and = U = r, yield the constraints 

1 = + Y?_ + £2 + Nl + N 2 _ , (4.7a) 

= iV x E_ - N- E x , (4.7b) 

and the evolution equations (|4.2|) . with Eq. (|4.3|) . and the decoupled evolution 
equation for Ei 1 , Eq. (j3. 13a|) . 

4.2.1. Taub invariant set. The Taub invariant set, a set of SH vacuum solutions 
of Bianchi Type-II, is defined as the subset 

= E x = iV x 

of Eqs. (j4.7|) . (|4.2j) and (j4.3j) . This leads to the reduced system of equations 

1 = E^. + Y?_ + Nl , (4.8) 



and 



C- X (l + £+) 9 t (l + S+) = (q - 2) (1 + S+) , (4.9a) 
C -1 (l + S+) = (g - 2) E_ - 2\/3 AT^ , (4.9b) 
C _1 (l + E+) 9tiV_ = (g + 2E+ + 2 V / 3S_) iV_ , (4.9c) 



with 



g = 2(E^ + El). (4.10) 
On account of employing Eq. ()4.8j) in Eq. ()4.1U|) . one finds from Eqs. ()4.9a() 
and (|4.9b|) that the Taub orbits are represented in the (E_|_E_)-plane of the 
Hubble-normalized state space by a 1-parameter family of straight lines, i.e., 

E_ = constant x (1 + E+) - v^ , (4.11) 

(see WE, Subsec. 6.3.2, for further details). 

4.2.2. Kasner invariant set. The Kasner invariant set, a set of SH vacuum 
solutions of Bianchi Type-I, is defined as the subset 

= N x = iV_ 

of Eqs. (|4.7|) . (|4.2j) and (|4.3|) . This leads to the reduced system of equations 

1 = E^ + E^ + E 2 X , (4.12) 

and 

C7- 1 (l + E + )a t (l + E+) = , (4.13a) 

C -1 (l + E+) a t E_ = 2v / 3S 2 , (4.13b) 

C- 1 (l + E+)a t E x = -2 V / 3S_E X , (4.13c) 

with 

9 = 2. (4.14) 

These equations define a flow on the unit sphere in E-space, with coordinates 
(E + , E_, E x ), which we refer to as the Kasner sphere. The Kasner orbits on 
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this sphere represent Kasner solutions with respect to a spatial frame that is not 
Fermi propagated along eo; such orbits are depicted in Fig. ^ The projections 
of the Kasner orbits into the (E + S_)-plane of the Hubble-normalized state 
space form a 1-parameter family of parallel straight lines, given by 

S + = constant . (4-15) 

The subset of Kasner orbits that represent Kasner solutions with respect to 
a spatial frame that is Fermi propagated along eo is given by [ cf . Eq. (|;-{.11J) ] 

= £x , (4.16) 

thus giving rise to a unit circle of equilibrium points in the (S+S_)-plane de- 
scribed by 

1 = E^- + , = <9 t X+ = d t ^- . (4.17) 
This set of equilibrium points forms the familiar Kasner circle, denoted by 

K. . 

The Kasner exponents in the equilibrium state are given by (WE, Subsec. 6.2.2) 

Pi = \0-- 2S +) . ^2,3 = I (l + S + ± v^X-) . (4.18) 

Note that according to an ordering given by pi < P2 < P3 and its permutations, 
K, subdivides into six equivalent sectors. The sectors are separated by special 
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points for which the variable pair (£+, £_) takes the fixed values 
Qi ■ (1,0) , T 3 






Q 3 

so that two of the Kasner exponents become equal. The points Q a correspond 
to Kasner solutions that are locally rotationally symmetric, while the points 
T a correspond to Taub's representation of the Minkowski solution (see WE, 
Subsec. 6.2.2, for further details). 

5. Asymptotic dynamics towards the singularity 

5.1. Description using metric variables. In this subsection we will briefly 
review the dynamical phenomenology of Gowdy vacuum spacetimes as described 
by the system of equations Q2.2[) for the metric variables A, P and Q; this was 
discussed earlier by Berger and Moncrief |Hj and Berger and Garfinkle [5] on the 
basis of extensive numerical experiments. In the following we reproduce some 
of their results, using LeVeque's clawpack 10 ; see Ref. j-ilj for background. 
Employing the "standing wave" initial data of Ref. [S], 

P(0,x)=0, P t (0,x) = 5cosx , 

Q(0, x) = cos x , Q,t(0, x) = , 

A(0,x) = , 

our integration proceeds on a 1-dimensional spatial grid representing the x- 
interval (0, 2ir), of step size 0.2 x 10~ 3 (thus there are 31416 spatial points; end 
points identified), with a time step size of 0.2 x 10~ 3 . 

In our analysis we are only concerned with the direction towards the singu- 
larity, which, with our choice of area time coordinate, corresponds to t — > +00. 
In this regime, because of the factors e~ 2t accompanying the spatial derivatives 
in Eqs. (|2.2|) . it is plausible that the evolution along different timelines (each 
of which being given by x = constant) decouples. As is well known, in the ap- 
proach to the singularity the metric variables P and Q (as functions of t and x) 
develop so-called "spiky features"; see Fig. for a plot at a fixed t. 

The apparent discontinuities in Q at certain values of x develop along excep- 
tional timelines for which Q t = 0. They are induced by local minima in P with 
P < 0. The latter are known to be coordinate (gauge) effects and are hence 
referred to by Rendall and Weaver [31] as "false spikes" . 

The sharp local maxima in P with P > 0, on the other hand, develop along 
exceptional timelines for which Q tX = 0. That is, Q exhibits smooth behavior 
along such timelines. The upward pointing peaks in P are known to have a 



ia 



This package is available from the URL: http : //www. amath. Washington. ed u/^claw/| 
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geometrical origin and are hence referred to by Rendall and Weaver [33] as 
"true spikes" . 

In terms of the hyperbolic space geometry, the false spikes correspond to 
points on a loop in H 2 , evolving under the Gowdy equations, which travel 
towards a point on the conformal boundary of H 2 chosen as infinity in passing 
to the model given by the metric 1)2 .5 J) . On the other hand, the true spikes 
correspond to points located on cusps of the loop, which travel at higher velocity 
than the rest of the loop. 

Note that while the coordinate derivative quantities and Pj are scale- 
invariant, P x , Qj and Q tX are not. Examples of geometrically defined quantities 
are the hyperbolic velocity and the hyperbolic norm of the spatial derivative, 
which, with u(t,x) representing the point in H 2 , are given by 

v(t,x) = \\d t u\\g H2 , e~* \\d x u\\ gii2 . 
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Since we are using a logarithmic time coordinate, these two quantities are scale- 
invariant by construction. One is inexorably led to consider the scale-invariant 
operators c\ and e~* d x . In terms of the dimensional time coordinate r = ^o e_t ; 
these operators are — t d T and rd x , respectively. 

The "energy density", A ; j, depicted at a fixed t in Fig. @J decreases more 
slowly at spike points. Since we are not using an adaptive solver, we are not 
able to numerically resolve the behavior of A f at spike points. With u(t, x) 
representing the point in H 2 , however, we have from the Gaufi constraint (|2..3a|) 

A, t (t, x) = v 2 (t, x) + e~ 2t \\d x u\ | 2 h2 , 

where the second term is expected to tend to zero pointwise. Therefore, we 
expect that A^ asymptotically tends to a nonzero limit pointwise, given by the 
asymptotic hyperbolic velocity which we define as 

v(x) := lim v(t,x) , (5.1) 

t— >+oo 

where the limit exists. As argued by Berger and Moncrief jSj and Berger and 
Garfinkle jS], the dynamics force the hyperbolic velocity to values v < 1 at late 
times, except along timelines where spikes form; see Fig. OD The reason max^ v 
appears to be less than 1 at late times in Fig. El is that the spikes have become 
so narrow that they are not resolved numerically. 

Rendall and Weaver [31] have shown how to construct Gowdy solutions with 
both false and true spikes from smooth solutions by a hyperbolic inversion 
followed by an Ernst transformation. In particular, they show that along those 
exceptional timelines, x = XtSpike; where (in numerical simulations) true spikes 
form as t — > +oo, the asymptotic hyperbolic velocity satisfies v(xtSpike) = 1 + s, 
with s £ (0, 1). On the other hand, along nearby, typical, timelines the limiting 
value is v = 1 — s. The speed at which v approaches v = 1 — s along the 
latter timelines decreases rapidly as x approaches a^tSpikc- ^ or ^ ne exceptional 
timelines, x = Xfs p ikc> where (in numerical simulations) false spikes form as 
t — ► +oo, the asymptotic hyperbolic velocity satisfies t>(xfs p ike) = s, again with 
s G (0, 1). Nearby, typical, timelines have limiting velocities v close to s. The 
regimes 1 < v{x t s P ike) < 2 and < £(xfs P ike) < 1 are referred to by Garfinkle 
and Weaver as the regimes of "low velocity spikes" . 

The true spikes with u(x t g p ike) £ (1)2) correspond, as has been known for a 
long time, to simple zeros of Q^ x . Higher velocity spikes were also constructed 
by Rendall and Weaver. However, these correspond to higher order zeros of 
Q jX , and are therefore nongeneric. 

5.2. Description using Hubble-normalized variables. Here we first de- 
scribe the behavior of the Hubble-normalized variables E+, £_, N x , S x , 
and U, when expressed in terms of the metric variables A, P and Q and their 
coordinate derivatives as in Eqs. (|3.18|) . under the numerical integration of the 
Gowdy equations (|2.2[) - (|2.3[) . We use the same "standing wave" initial data as 
in Subsec. 15.11 The spatial variation of the values of the Hubble-normalized 
variables at a fixed t is depicted in Figs. Ifi hlll Note, in particular, the scales 
along the vertical axes in these figures. 
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Fig. 6. £+ at t = 10. 
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Fig. 7. U at t = 10. 
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Fig. 8. XL at t = 10. 
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Fig. 9. iV x at t = 10. 
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Fig. 10. E x at t = 10. 



Fig. 11. N- at t = 10. 



As we see in Figs. |T| and U and -/V x become uniformly small as t increases. 
On the other hand, X 2 tends to the value 1 almost everywhere except along 
those timelines where either a false or a true spike is being formed; see Fig. 1121 
Further, while X x and iV_ do not necessarily become small individually as t 
increases, see Figs. ITUl and 1111 their mutual product A^_S X , nevertheless, does; 
see Fig. 1131 Relating these observations and the analytic result (derivable from 
Subsec. 13. If) that 

l^i 1 ! < 4e _t , (5.2) 
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so E\ — ► uniformly as i — ► +00, to the system of equations (|3.12|) - (|3.14|) 
suggests that, in the approach to the singularity, for Gowdy vacuum space- 
times the ultra-strong spacetime curvature phenomenon of asymptotic silence, 
described in Refs. [2U1 §5.3] and jSHl §4.1], sets in. Asymptotic silence refers 
to the collapse of the local null cones onto the timelines as t -> +00, with 
the physical consequence that the propagation of gravitational waves between 
neighboring timelines is gradually frozen in. Based on these considerations, it is 
natural to expect that near the singularity the spatial derivatives of the Hubble- 
normalized connection variables become dynamically irrelevant, and that thus, 
in a neighborhood of the silent boundary of the Hubble-normalized state space, 
the dynamics of Gowdy vacuum spacetimes is well approximated by the dy- 
namics on the silent boundary. More precisely, it is well approximated by the 
system of equations (|3.12bj) . (|3.12cjl . lJO|> . and (flTTaj) . but note that the 
Hubble-normalized variables are now x-dependent. This aspect of the asymp- 
totic dynamics may also be viewed as a consequence of the friction term r" 1 d T 
in Eq. (|2.6j) . In Figs. E] an dE3 w e show the projections into the (S + S„)-plane 
of the Hubble-normalized state space of orbits for (i) the evolution system on 
the silent boundary, Eqs. (|4.2j) . and (ii) a typical timeline of the full Gowdy 
evolution equations (|3.13j) . It is apparent that qualitative features of the full 
Gowdy evolution system along an individual timeline are correctly reproduced 
by the evolution system on the silent boundary. 

To gain further insight into the dynamics towards the singularity, we will now 
investigate the local stability of the Kasner circle IC on the silent boundary of 
the Hubble-normalized state space for Gowdy vacuum spacetimes, thus probing 
the closest parts of the interior of the state space. 

5.2.1. Linearized dynamics at the Kasner circle on the silent boundary. We 
recall that on the silent boundary the Hubble-normalized variables are x-depen- 
dent. Consequently, the Kasner circle IC on the silent boundary is represented 
by the conditions 

= Ex 1 = S x = iV x = N- = U = r , £+ = £+, £_ = £_, (5.3a) 

1 = ±l + t 2 _ , q = 2 , (5.3b) 
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Fig. 14. Projection 
into the (£+£_)- 
plane of an orbit de- 
termined by the SB 
system (|4~2|) . 



Fig. 15. Projection 
into the (£+£_)- 
plane of a Gowdy or- 
bit along the typical 
timeline x = 2.7. 



where we introduce the convention that "hatted" variables are functions of the 
spatial coordinate x only. Thus, linearizing the Hubble- normalized Gowdy evo- 
lution equations ()3. 13|> on the silent boundary at an equilibrium point (S+, S_) 
of fC yields (with C = — \) 



d t N x 
dt^x 



Ei 1 , 



i 



2 + 



:i + s_ 



i + 



(l + E- 



(5.4a) 
(5.4b) 

(5.4c) 
(5.4d) 
(5.4e) 
(5.4f) 



here and denote the deviations in S + and £_ from their equilibrium 
values and £_ on /C. Let us define an x-dependent function 



(5.5) 
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corresponding to a particular function used in Refs. j^Hj an d |34| . Then, inte- 
grating Eqs. 1)5. 4 J) with respect to t, we find for E\ , and iV x 



= 0(e - *) , (5.6a) 

= 0(e - ') , (5.6b) 

= 0(e - *) , (5.6c) 

N x = N x e _< + 0(ie _ *) . (5.6d) 



That is, as i — ► +oo, the solutions for E\\ 5T, + and <5£_ decay exponentially 
fast as e - *, while iV x decays exponentially fast as ie - *, independent of the 
values of (£+,£_) [and so k(x)] along an individual timeline. Consequently, 
these variables are stable on KL. The solutions for the variable pair (E x ,iV_) 
associated with the "x -polarization state", on the other hand, are given by 

S x =S x e- fc ^ , (5.7a) 

jV_ =iV_ e -[i-« , (5.7b) 

(see also Eqs. (140) and (141) in Ref. EDI)- 11 Thus, 

(1) S x is an unstable variable on K, as t — ► +oo whenever fc(rr) < 0; this 
happens along timelines for which the values of £_) define a point 
on the upper semicircle of /C (unless, of course, = E x along such a 
timeline), and 

(2) iV_ is an unstable variable on KL as t — > +oo whenever /c(x) > 1; this 
happens along timelines for which the values of S_) define a point 
either on the arc(TiQ3) C /C or on the arc(QsT2) C K. (unless, of course, 
= iV_ along such a timeline). 

This implies that in the approach to the singularity, the Kasner equilibrium 
set /C in the Hubble-normalized state space of every individual timeline has a 
1- dimensional unstable manifold, except when the values of (S+,2_) coincide 
with those of the special points Q a and T a , or correspond to points on the 
arc(T2<5i) C K, which is stable. 12 

Close to the singularity (and so to the silent boundary of the state space), 
it follows that along those timelines where (S + ,E_) take values in the Ex- 
unstable sectors of /C, a frame transition orbit of Bianchi Type-I according to 
Eq. 1)4. 15J) to the S x -stable part of K, is induced (corresponding to a rotation 
of the spatial frame by ir/2 about the ei-axis). Notice that timelines along 
which = S x do not participate in these frame transitions and so are "left 
behind" by the asymptotic dynamics. This provides a simple mechanism for 
the gradual formation of false spikes on T 3 . As, by Eq. ()3.11j) . E x quantifies 
the Hubble-normalized angular velocity at which e.<i and e% rotate about e.\ 
(the latter being a spatial frame gauge variable), false spikes clearly are a gauge 
effect. 

^Altogether, the linear solutions contain four arbitrary 27r-periodic real-valued functions 
of x, namely k, N x , E x and AT_. 

1 9 

For the polarized invariant set, on the other hand, here included as the special case 
= E x = N-, the entire Kasner equilibrium set K, is stable. 



21 



L. ANDERSSON, H. VAN ELST, AND C. UGGLA 



On the other hand, along those timelines where (£+,£_) take values in 
the iV_-unstable sectors of IC, a curvature transition orbit of Bianchi Type-II 
according to Eq. 1)4.11(1 to the iV~_-stable part of IC is induced; this is shown in 
Fig. 1161 Notice that timelines along which = iV_ do not participate in these 
curvature transitions and so are "left behind" by the asymptotic dynamics. 
This provides a simple mechanism for the gradual formation of true spikes 
on T 3 . As iV_ presently relates to the intrinsic curvature of T 3 , true spikes are 
a geometrical effect. 

It was pointed out in Ref. |38l §4.1] that it is natural to expect from the 
consideration of asymptotic silence that 

lim E 1 1 d x Y = (5.8) 

t— >+oo 

is satisfied along typical timelines; in our case Y = (£+, S_, -/V x , £ x , iV_) T . 
This limit is certainly attained for each of £+, £_ and A^ x when we use the 
linear solutions of Eqs. (|5.6|) for these variables which are valid in a neighbor- 
hood of IC on the silent boundary. For the "x -polarization variables" (£ x , N-) 
an analogous substitution from Eqs. (|5.7|) leads to 

^S x oc (d x ± x - t± x d x k{x)) e-^Wl* , (5.9a) 
E x x d x N- oc (d x N- +tN^ d x k(x)\ e~ [2 - k{x)]t . (5.9b) 

Here a number of cases arise when we consider the behavior of these spatial 
derivative expressions in the limit t — ► +oo. Firstly, both decay exponentially 
fast along those timelines for which the dynamics has entered the regime < 
k(x) < 1. As discussed in some detail below, this happens asymptotically along 
typical timelines. Secondly, along timelines where false spikes (0 = S x ) form, 
both decay exponentially fast when the dynamics is confined to the regime 
— 1 < k(x) < 0. Thirdly, along timelines where true spikes (0 = iV_) form, 
both decay exponentially fast when the dynamics is confined to the regime 
1 < k(x) < 2. Outside the regime — 1 < k(x) < 2, one or the other of the 
spatial derivative expressions (|5.9|) may grow temporarily along timelines, until 
one of the unstable modes (|5.7J) drives the dynamics towards a neighborhood 
of a different sector of IC. Also, outside the regime — 1 < k(x) < 2, spikes of 
higher order may form along exceptional timelines, through choice of special 
initial conditions. 

In recent numerical work, Garfinkle and Weaver |2J investigated spikes with 
k{x) initially outside the regime — 1 < k{x) < 2 and reported that they generally 
disappear in the process of evolution. It has been suggested by Lim [35] that so- 
called spike transition orbits could provide an explanation for this phenomenon. 
Unfortunately, no analytic approximations are available to date and further 
work is needed. Nevertheless, numerical experiments so far indicate that the 
limit (|5.8|) does hold for timelines with nonspecial initial conditions |32j . 

According to a conjecture by Uggla et al jSHJ §4.1], the exceptional dynamical 
behavior towards the singularity along those timelines where either false or 
true spikes form does not disturb the overall asymptotically silent nature of the 
dynamics, as in the (here) asymptotic limit t — ► +oo any spatial inhomogeneity 
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associated with the formation of spikes is pushed beyond the particle horizon 
of every individual timeline, and so cannot be perceived by observers traveling 
along the timelines. 13 

The present consideration provides a natural basis for classifying spikes as 
"low velocity" or "high velocity", according to Rendall and Weaver [31] and 
Garfinkle and Weaver [2T]. By Eqs. (f53|) . 

(1) "low velocity" arises (a) for false spikes when = S x and —1 < k(x) < 
0, (b) for true spikes when = iV_ and 1 < k{x) < 2, 

(2) "high velocity" arises (a) for false spikes when = £ x = d x T, x and 
k(x) < — 1, (b) for true spikes when = iV_ = d x N_ and k(x) > 2. 

Due to the strong constraints that need to be satisfied, (persistent) high velocity 
spikes are nongeneric. 

For fixed x, the definition of the function k(x) given in Eq. (|5.5|) provides 
an alternative parametrization of the Kasner circle. This is of some interest 
since, according to Kichenassamy and Rendall |29| p. 1341], the asymptotic 
hyperbolic velocity v(x) defined in Eq. (|5.1j) and k{x) are related by 

v(x) = \k(x)\ . 

In this parametrization of /C, the special points Q a and T a correspond to the 
values 

Qi : k = , T 3 : k = - 1 , 

Q2 : k = — 3 , T\ : k = 00 , 

Q 3 : k = 3 , T 2 :k=l . 

In addition, the value k = 2 corresponds to = (— 1/7, — 4\/3/7). The 

marked points on the Kasner circles in Figs. El an dEl are (from the left) the 
points with k = 00, 2, 1, 0. 

Returning to the discussion of the asymptotic dynamics as t — > +00, as briefly 
mentioned above our numerical simulations suggest that the product between 
the "x -polarization variables" (S x , N-), which span the unstable directions on 
the unstable sectors of the Kasner circle (and are also responsible for the spike 
formation process), satisfies the limit 

lim iV_£ x = . (5.10) 

This behavior, a special case of Eq. (101) in Ref. |3H]) is further reflected in the 
linear solutions given in Eqs. (|5.7j) . which, for t — > +00, yield exponentially fast 
decay of -/V_£ x as e~ l . In the same regime, A^ X S_ decays exponentially fast as 
te~ f , so that, by Eq. ()3.12c|) . U decays as ie -4 , and so does r by Eq. (|3.12a|) . 
The linear analysis thus strengthens the emerging picture of the dynamical 
behavior along typical timelines of Gowdy vacuum spacetimes in the approach 
to the singularity. 

13 In unpublished work, Woei Chet Lim and John Wainwright have constructed an explicit 
exact solution with a true spike for which the conjecture has been shown to be analytically true. 
Moreover, unpublished numerical work by Lim suggests that for Gowdy vacuum spacetimes 
this is true in general. Thus, there exists evidence for asymptotic silence to hold even when 
spatial derivatives blow up due to the formation of spikes. 
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Fig. 16. Curvature 
transition orbit of 
Bianchi Type-II. 



Fig. 17. Snapshot 
of the solution at 
t = 12, projected 
into the (£+£_)- 
plane. Frame tran- 
sition and curvature 
transition orbits are 
indicated. 



Asymptotic silence holds for Gowdy vacuum spacetimes since Ei 1 — ► uni- 
formly as t — > +oo; see Eq. Q5.2JI . Our numerical calculations indicate that 
iV_X x — > and iV x X„ — ► uniformly (see Fig. I13j) so that U — ► and r — > 
uniformly as t — > +oo. It follows that the dynamics of Gowdy vacuum space- 
times are increasingly well approximated by the dynamics on the silent bound- 
ary. The latter lead to a characteristic interplay between S x -induced frame 
transitions and iV_-induced curvature transitions on the Kasner circle associ- 
ated with each timeline; they give rise to a finite sequence of transitions from 
one unstable position on K, to another one, until, eventually, a position on the 
stable arc(T2<5i) of IC is reached. This behavior is shown clearly in Figs. llBh 
1171 The interplay between the frame transitions and the curvature transitions 
provides a simple mechanism for forcing the function k(x) to satisfy 

\k(x)\ < 1 , 

along typical timelines. Understanding the reduction of the value of the hyper- 
bolic velocity v(t,x) to the regime v(x) = \k(x)\ < 1 along typical timelines in 
the evolution of Gowdy vacuum spacetimes towards the singularity is stated as 
an open problem by Kichenassamy and Rendall |291 p. 1341]. 



5.3. Behavior of Weyl curvature. We now evaluate the asymptotic behavior 
of the Hubble-normalized electric and magnetic Weyl curvature (relative to eo) 
for Gowdy vacuum spacetimes in the so-called "low velocity" regime < k(x) < 
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1, i.e., for typical timelines with the dynamics in a small neighborhood of the 
stable arc(T2<5i) C K,. Substituting the linear solutions 1)5.6)) and 1)5.7(1 into the 
Weyl curvature formulae given in Add. |A~T1 we obtain 



£+ 


6\l-k 2 (x)} , nl f , , nl 


-2[i-k(x)]u + , e -2k(x)u ^ 


(5.11a) 


■H+ 


= 0( e - [1 - fc(x)]i ) + 0(te~ [1+fc(x)lt ) 




(5.11b) 


£_ 


2 V '3fc(x)[l-fc 2 ( a ;)] _ t 
[3 + fc 2 (x)] 2 + 1 J 


+ ( e -2[l-*(*)]t) ; 


(5.11c) 


n x 






(5. lid) 


s x 


= 0{e~ k{x)t ) , 




(5.11e) 


H- 


= 0(e~ [1 - k{x)]t ) + 0(te- [l+k{x)]t ) 




(5.11f) 



In the limit t — > +oo, the magnetic Weyl curvature decays exponentially fast, 
so that for the electric Weyl curvature the Kasner limit is attained pointwise. 

6. Concluding remarks 

In Go cosmology, as in the Gowdy case, the evolution equations on the silent 
boundary for the Hubble-normalized dimensionless coordinate scalars (i.e., all 
variables except the frame variables E a l ), given in Ref. [381 §2], are identical 
to those for the same quantities when describing SSS and SH models. The 
reason for this is the following. Dimensionless scalars of SSS and SH models 
are by definition purely time-dependent (in symmetry adapted coordinates), 
and thus the spatial frame derivatives that appear in the equations for the 
Hubble-normalized dimensionless coordinate scalars drop out, which leads to 
a decoupling of the frame variables E a l . The resulting system of equations is 
clearly the same as that obtained by setting the spatial frame derivatives to 
zero, as done on the silent boundary. 

In the identification of the equations on the silent boundary with the SSS 
and SH equations, the normalization is essential; it is only when dealing with 
scale-invariant dimensionless variables that one obtains a direct correspondence. 
As discussed by Eardley [T5| . a SSS geometry is related to a SH geometry by 
a spatially dependent conformal transformation; by introducing dimensionless 
variables the conformal factor drops out and one obtains SH quantities. The 
silent boundary of the Go state space consists of a union of the most general SSS 
and SH models, notably SH Bianchi Type- VIII and Type-IX models and self- 
similar models of class D (see Ref. |151 §3] for notation). 14 To obtain a detailed 
proof of the SSS/SH silent boundary correspondence, one can start with the 
SSS and SH symmetry adapted metric representation given by Eardley ^3 §3], 

14 The joint SSS and SH character of the silent boundary is in accordance with the attempt 
in Ref. [38] to heuristically explain the dynamical importance of the silent boundary: ultra- 
strong gravitational fields associated with typical spacelike singularities collapse the local null 
cones, which thus prevents information propagation yielding local dynamics. SSS spacetimes 
are conformally SH, and thus also SSS dynamics are purely local. Moreover, the "pathologi- 
cal" spatial properties of SSS spacetimes are irrelevant; only their temporal properties are of 
interest in the approach to a spacelike singularity. 
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and then derive the correspondence through a direct comparison between the 
SSS/SH equations and those on the silent boundary. 

Gowdy vacuum spacetimes form an invariant set of the general Go state 
space. The intersection of this invariant set with the SSS and SH subsets in 
turn yields an invariant subset which is described by the equations on the silent 
boundary of the Gowdy state space. 

6.1. AVTD and BKL from the point of view of the silent boundary. 

Discussions of AVTD behavior near spacetime singularities and of the BKL 
conjecture on oscillatory behavior of vacuum dominated singularities have of- 
ten been associated with claims that the asymptotic dynamical behavior in 
the direction of the singularity of spatially inhomogenous spacetimes is locally 
like that of SH models. Indeed, BKL take as a starting point for their anal- 
ysis explicit SH solutions. They replace the integration constants by spatially 
dependent functions and perform a perturbation analysis in the asymptotic 
regime around the resulting ansatz. The dynamical systems approach offers 
some justification for this procedure. 15 

Due to the SSS/SH silent boundary correspondence, the reduced system of 
equations on the silent boundary leads to the same solutions for the Hubble- 
normalized coordinate scalars in both cases. However, in contrast to the true 
SSS/SH cases, the constants of integration become spatially dependent func- 
tions on the silent boundary. Then, in the true SSS/SH cases, the metric is ob- 
tained in a subsequent step, by integrating the decoupled system of equations 
for the frame variables E a % . However, the silent boundary (where = E a l ) 
is unphysical, since it corresponds to a degenerate metric. In order to obtain 
an expression for the metric which is valid in the asymptotic regime, one is 
forced to go into the physical interior part of the state space where the E a l are 
nonzero. By perturbing around = E a l , starting from a seed solution to the 
reduced system of equations on the silent boundary, approximate solutions to 
the full system can be constructed. To lowest order, the asymptotic approxi- 
mation so obtained is identical to the corresponding true SSS/SH solution, but 
with integration constants replaced by spatially dependent functions which in 
turn are restricted by the Codacci constraint. 

In terms of the discussion in this paper, the AVTD leading order solution 
for Gowdy vacuum spacetimes is produced by the above procedure using a seed 
solution taking values in the Kasner subset on the silent boundary (note that 
VTD is associated with the entire silent Kasner subset, while AVTD is associ- 
ated only with the stable part(s) of the Kasner circle). Indeed, the linearization 
at the stable arc of the Kasner circle yields the starting point for the analysis 
of Kichenassamy and Rendall Similarly, BKL in their analysis consider 

seed data in the Taub subset (vacuum Bianchi Type-II) on the silent bound- 
ary, in addition to the Kasner subset. Thus, the AVTD and BKL analysis has 

15 An alternative point of view on the oscillatory approach to the singularity is provided 
by the method of consistent potentials (see Ref. Q f° r a discussion of the T 2 -symmetric 
case). Further, Thibault Damour, Marc Henneaux and collaborators have recently developed 
a picture of the asymptotic behavior of stringy gravity in D = 10 spacetime dimensions, which 
leads to consideration of hyperbolic billiards (see Refs. 1 1 31 and |14| . and references therein). 
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considered only part of the past attractor conjectured in Ref. while on the 
contrary, the results in the present paper indicate that it is essential to consider 
the entire silent boundary in order to understand the approach towards the 
singularity. 

6.2. Gowdy case. We conclude by discussing the picture that emerges in the 
Gowdy case. In view of the non-oscillatory nature of the system of equations 
on the silent boundary, we expect Gowdy vacuum spacetimes to have a non- 
oscillatory singularity in the sense that lim^ +00 X(t, x) = X(x) exists, and 
that X takes values only on the Kasner circle fC on the silent boundary. More- 
over, the local stability of the arc(T2(5i) C K. suggests that it is a local attractor, 
i.e., 

A = arc(T 2 Qi) C K = A^ owdy . 

This was indeed proved by Ringstrdm |36l I37j and Chae and Chrusciel [5] : if 
X(t,x) is sufficiently close to A, in a suitable sense, then lim^+oo X(t, x) = 
X(x), with X(x): S 1 — > A a continuous map, and the limit is in the uni- 
form topology. Chae and Chrusciel [H] have also proved that for any Gowdy 
vacuum spacetime with smooth initial data, there is an open and dense sub- 
set O C S l , such that for x € O, lim f ^ +00 X(t, x) = X(x), X(x) G A, and 
X is smooth at x. Furthermore, they showed that, for any closed F C S 1 
with empty interior, a solution could be constructed with (false) spikes on F. 
By applying Gowdy-to-Ernst transformations, these can presumably be turned 
into true spikes with high velocity. This shows that the detailed asymptotic 
behavior of Gowdy vacuum spacetimes can be quite complicated. A very useful 
further step in the mathematical analysis of Gowdy vacuum spacetimes would 
be to prove that it is always the case that lim i ^ +00 (S_A r x ) 2 + (X x iV_) 2 = 0. 
A limit of this kind reflects that in the asymptotic regime the propagation of 
gravitational waves (represented by the "+ -polarization variables" (£_,iVx) 
and the "x -polarization variables" (E x ,iV_), respectively) becomes dynami- 
cally insignificant. 

If one considers generic smooth initial data, the picture should simplify con- 
siderably. For example, since a generic smooth function has at most a finite 
number of zeros on a finite interval, and true spikes correspond to zeros of 
Q tX , one expects that generic solutions have at most a finite number of true 
spikes. An analogous argument indicates that a generic solution has at most a 
finite number of (true or false) spikes. High velocity spikes are expected to be 
nongeneric since they correspond to higher order zeros of Q tX . Therefore, one 
expects that a generic solution will have a finite number of true spikes, all with 
velocity in the interval (1,2). 

Due to the existence of spikes for general Gowdy vacuum spacetimes, A can- 
not be a global attractor. For a Gowdy solution with spikes, 

lim sup d(X(t,x),JC) ^ , 

*-H-°° x6 Si 

even though we expect that the pointwise limit is a map to 1C. This makes 
it clear that the notions of "asymptotic state" and "attractor" depend, in the 
spatially inhomogenous case, on the choice of topology. 



30 



L. ANDERSSON, H. VAN ELST, AND C. UGGLA 



The numerical work in this paper indicates that the variety 

U = {S_iVx = E x iV_ = 0} 
is the uniform attractor for Gowdy vacuum spacetimes in the sense that 

lim sup d(X(t,x),U) = . 

As seen from the discussion in Subsec. 15.21 the silent boundary dynamics on U 
explain the major evolutionary features along typical Gowdy timelines. 

The state space framework for spatially inhomogeneous cosmology based on 
the dynamical systems approach offers a common ground where many scattered 
results and conjectures may be collected and put into a broader context. For 
example, earlier claims and results as regards SH models can now be regarded 
as results determined by the structure of the silent boundary of the Hubble- 
normalized state space. This pertains, e.g., to work by BKL, conjectures and 
proofs in WE, and proofs by Ringstrom. The latter showed in Ref. [351 that 
the Bianchi Type-II variety is a past attractor for nontilted SH perfect fluid 
models of Bianchi Type-IX. 16 

Finally, earlier results, and the results in this paper, suggest that the silent 
boundary, perturbations thereof, and the formation of spikes, in a Gi as well as 
in a Go context, are worthy areas of exploration, and are likely to offer exciting 
hunting grounds in the coming years. 
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Appendix A. Hubble-normalized Weyl curvature 

A.l. Variables. The Hubble-normalized electric and magnetic Weyl curva- 
tures (relative to eo) for Gowdy vacuum spacetimes can be easily obtained by 
specializing the relations given in Ref. ^Hj. The "+", "— " and "x" variables 



This attractor organizes the past asymptotic dynamics and provides an explanation for 
the Kasner-billiard-like behavior of nontilted SH models of Bianchi Type-IX in the approach 
to the singularity. The same variety is expected to be the past attractor for nontilted SH 
models of Bianchi Type- VIII. A candidate attractor variety has been identified for nontilted 
SH perfect fluid models of Bianchi Type- VII x / 9 by Hewitt el al 1251. This is the only nontilted 
SH model in class B that exhibits oscillatory dynamical behavior into the past. 
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are defined analogous to the shear rate variables in Eqs. Q2.13JI . 

5+ = | (N x + Nl ) + i (1 + E + ) E+ - 1 (S 2 „ + S 2 X ) , (A.la) 

H+ = -A_£_ - iV x S x , (A.lb) 

£_ = i (i^i 1 ^ - r) iV x + A N 2 _ + i (1 - 2£+) , (A.lc) 

W x = i (^i 1 «9 X - r) £_ - A ]\r_E x - iV x S+ , (A. Id) 

f x = - i (^i 1 9 X - r) A_ + Ar x AT„ + I (i _ 2S+) E x , (A.le) 

W- = - \ (Ei 1 d x -r)E x -^= iV_E_ - iV_E+ . (A.lf) 

A. 2. Spatial scalars. 

5 a/3 5 a/3 = 6(5 1 + 5! + 5 X ) , (A.2a) 

n aP n aP = 6(H% + ni + n 2 x ) , (A.2b) 

5 Q /3W Q/3 = 6(5+ W+ + 5_W_ + 5 X W x ) , (A.2c) 

SjSfPS^ = - 65+ (5+ - 35! - 35 2 ) , (A.2d) 

HcfHpm^ = - m+ {n 2 + - m 2 _ - m\ ) , (A.2e) 



£ a p Hp f 'H 1 a = - 6 [5+ (H% - H 2 _ - Til) - - 2£ X H + H X ] , 

(A.2f) 

njEfPe^ = -6[n+ (5 2 - el -e x )- 2w_5+5_ - 2n x e + e x ] . (A.2 g ) 

A. 3. Spacetime scalars. 

Ji = 8(e a pe al3 - H aP H aP ) , (A.3a) 

J 2 = - 165 a/3 W a/3 , (A.3b) 

T 3 = - 16(5 Q %T5 7 Q _ zsjHpm^) , (A.3c) 

j 4 = i6(njn^n 7 a - mje^ej 1 ) . (A.3d) 

2+ is the Hubble- normalized Kretschmann scalar. 
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